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g ■ Abstract 

^C) . We investigate the muon anomalous magnetic dipole moment (MDM), 

| the muon electric dipole moment (EDM) and the lepton-flavour- violating 

■ decays of the t— lepton, t — * /17 and t — * 3/i, in the CP-violating Min- 

imal Supersymmetric Standard Model (MSSM) with nonuniversal soft- 
supersymmetry breaking. We evaluate numerically the muon EDM and 
the branching ratios B(r — ► fyf) and B{t — ► 3/i), after taking into ac- 
count the experimental constraints from the electron EDM and muon 
MDM. Upon imposition of the experimental limits on our theoretical 

' predictions for the aforementioned branching ratios and the muon MDM, 

jj^ ■ we obtain an upper bound of about 1CP 23 e-cm on the muon EDM which 

lies well within the explorable reach of the proposed experiment at BNL. 



PACS numbers: ll.30.Er, 12.60.Jv, 23.40.Bw 



1 Introduction 



Recently, the Muon— — 2 collaboration has reported the world average experimental value on the muon 
anomalous magnetic dipole moment (MDM) [1]: 

= \{9n - 2 ) = (H659203 ± 8) x 10~ 10 , (0.7 ppm). 
In the framework of the standard model (SM), the value of a^(SM) is currently evaluated to be [2] 

a^(SM) = (11659177 ± 7) x 10" 10 , (0.6 ppm). 
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The experimental value differs from the SM prediction by 1.6 standard-deviation (1.6a). Even though the 
1.6cr deviation is not very serious, this gap might be filled up by a contribution from new physics beyond the 
SM. It seems that a weak-scale new physics would fix the discrepancy [3, 4]. In the framework of the SM, the 
contribution to a M is traditionally divided into several pieces 



QED , hadronic . EW 



The QED loop effects have already been computed to high orders [5, 6]. A thorough analysis on hadronic 
contributions to the muon anomalous magnetic dipole moment is presented in Ref.[7]. At the one- loop level, 
the contribution of the standard model is formulated as [8, 9, 10, 11, 12] 



ew 5 Gpm , 
a,, 



M 3 Sv^vr 2 



1 , .9,9 „ , m2 , 



1 + ^-K)+°(: 



5 v- -W ' - v m 2 • 
w 



The two- loop EW-sector contributions to are also discussed in Ref.[13]. Provided the recent measurement 
of a M is taken to be a signal of new physics beyond the SM, whose corrections to are extensively discussed 
in literature, some authors [14] have analyzed the muon anomalous magnetic dipole moment in the MSSM. 
Considering the possible CP- violation phases, Ibrahim et.al made a similar analysis in the N = 1 supergravity 
model[15]. A systematic analysis on lepton-flavor-violating processes and a^— value within the framework of 
supersymmetry see-saw mechanism was given by Hisano ei.aZ.[16]. Involving the coupling of the second- 
generation with the third generation superparticles, the contributions of R-parity violation to have been 
evaluated [17] and discussions on in the supersymmetric GUTs are also made by some theorists [18]. A 
comparative study on in various supersymmetric models has also been presented in some recent works [19]. 
Provided one of the neutral Higgs bosons is of small mass, the authors of Ref.[20] calculated the anomalous 
magnetic dipole moment of muon in the two-Higgs doublet model. An analysis of the muon anomalous 
magnetic dipole moment in other extensions of the standard model has been given [21]. Alternatively, we will 
apply the effective Lagrangian [22, 23] to analyze the muon anomalous magnetic dipole moment in the CP- 
violating MSSM with nonuniversal soft-supersymmetry breaking, i.e., in this interaction which violates the 
lepton flavor conservation is mediated by the nonuniversal soft-supersymmetry breaking parameters. It is well 
known that the lepton-flavor-violating decays are also ideal processes to detect possible 'new' physics beyond 
the SM. So far, the experiments have not found any substantive evidence of such processes yet, instead, the 
experimental observation only sets upper bounds on those decay branching ratios, for example B(t — > /ry) 
and B(t —>■ 3//) [26]. Obviously, any new physics must be constrained by these bounds. 

In this work, we investigate the lepton-flavor-violating decays r — > ^7, r — ► 3/x, and the muon MDM, 
EDM in the framework of the MSSM with the nonuniversal soft-supersymmetry breaking. In the super- 
symmetric theories, there are many new physical CP-violating phases that are absent in the SM. Con- 
sidering the renormalization condition on CP-odd Higgs[27], we can choose the //—parameter 1 in the su- 
perpotential, and set the non-diagonal elements of the bilinear soft-supersymmetry breaking parameters 
m ^/7' m Q//' m u/7' m r)"' ^ ^ ^ anc ^ so ^ trurnear couplings A u , A d , A 1 with the physical CP 
phases after properly redefining the fields in the theory 2 . Up to one-loop order, those CP- violating phases 
induce the mixing among the CP-even and CP-odd Higgs [28, 29] and modify the Higgs boson couplings to 
the up- and down- quarks, and to the gauge bosons drastically [29]. The current experimental lower bound 
on the mass of the lightest neutral Higgs can be reduced to 60 GeV [29]. At present, the experimental upper 
bound on the electron EDM is set [26]: |d e | < 0.5 x 10~ 26 e • cm. In order to rationally predict the muon 

1 Please be noted, here we use the fx— parameter following the literature, hope that it would not cause any confusion with the 
muon which is sometimes written as fx. 

2 At the Lagrangian level, all the couplings may be complex. However some phases are un-physical, for evaluating the physical 
processes, they are not necessary and we can remove those phases by redefining the wavefunction as — » e^ty. This step is the 
same as to define the physical CP phase of the CKM matrix elements in the SM case. 
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EDM, we need to take the electron EDM as a rigorous constraint into account. It is well known that the 
two-loop Barr-Zee-type diagrams [24] may also give a large contribution to the electron EDM [25], thus in our 
discussion, we include the relevant two-loop Barr-Zee-type contributions to the EDM of charged leptons . 

Here, we adopt the notation of Ref.[30], the relevant nonuniversal soft-supersymmetry breaking terms and 
Feynman rules can also be found in Ref.[30]. Our paper is organized as follows. In Sec. II, we introduce 
the CP-violating MSSM with nonuniversal soft-supersymmetry breaking. In Sec. Ill, we analyze the loop- 
corrections to the e J e / 7 effective vertex. The muon anomalous magnetic dipole moment and the decay width 
of r — > hj in the supersymmetric models are eventually formulated. The r — ► 3/x is analyzed in Sec. IV. 
Within the experimentally allowed range for the concerned parameters, our numerical analysis is presented 
in Sec.V. Upon imposition of the experimental limits on the theoretical predictions of the branching ratios 
and the electron EDM and muon MDM, we obtain an upper bound on the muon EDM. Then we will make a 
brief summary about the method and model we employ in this work and discuss the obtained results in the 
last section. The tedious formulae are collected in appendices. 



2 The MSSM with nonuniversal soft-supersymmetry breaking 

The most general form of the superpotential which has the gauge invariance and retains all the conservation 
laws of the SM is written as 

W = fietjHlHj + eijti^Hmk 3 + e ij h d IJ HlQ I j D J + e^H 2 ^ . (1) 

Here H 1 , H 2 are the Higgs superfields; Q 1 and L 1 are quark and lepton superfields in doublets of the weak 
SU(2) group, where 1=1, 2, 3 are the indices of generations; the rest superfields U 1 , D 1 and R 1 are the quark 
superfields of u- and d-types and charged leptons in singlets of the weak SU(2) respectively. Indices i, j are 
contracted for the SU(2) group, and h!, h u,d are the Yukawa couplings. To break the supersymmetry, the 
nonuniversal soft-supersymmetry breaking terms are introduced as 

Coft = -™ 2 Hl HtH} - m 2 H2 H?*H? - m^^lf - m 2 RlJ R I *R J - m 2 QiJ Q\*Qi ~ m^U 1 *^ 

-m 2 Dij b u D J + (miA B Ai + m 2 \ i A \ i A + m 3 X a G X a G + h.c.) + [^BeijHjH 2 + e^^ZJiF 
+e ij A d j H}Q I j D J + eijA^HfQfi' + h.c] , (2) 
where m 2 , , m 2 „ , m 2 , m 2 , m 2 , m 2 and m 2 are the square masses of the superparticles, 

H 1 H 2 LU R.U Q 1 - 7 U ij D j J 

m3, m\ denote the masses of X G (a = 1, 2, ••• 8), A^ (i = 1, 2, 3) and A#, which are the 

SU(3) x SU(2) x U(l) gauginos. B is a free parameter in unit of mass. In jC so ft, the nonuniversal terms are: 
(a): in the bilinear couplings: m 2 , , m 2 T , m 2 , with I ^ J; whereas for I=J, m 2 r . m 2 T . m 2 rr are the uni- 

v ' L IJ U IJ D IJ L n U 11 D n 

versal soft terms; (b): in the trilinear couplings: A 1 ^, A u iV A d with I ^ J are the nonuniversal part, whereas 
as I=J, A l n , A™ z , A d n are the universal part. A l i7 , AJ j: A^ (I, J = 1, 2, 3) are the soft-supersymmetry 
breaking parameters that result in mass splitting between leptons, quarks and their supersymmetric partners. 
With the soft-supersymmetry breaking terms in Eq.(2), we can study the phenomenology in the minimal 
supersymmetric extension of the standard model (MSSM). The resultant 6x6 square-mass matrix of the 
charged scalar leptons is written as 



>r , m 2 Lij + (m 2 eI + (± + s 2 w )co S 2Pm 2 )5 IJ -m c i M * tan (35 IJ + 2 -^LcpA l u 

" >L "I -m e i tan J + ^fw C/3 ^* m 2 RlJ + (m^ + s% cos 2/?m|)5 /J ! ' 
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while the 3x3 sneutrino square-mass matrix is expressed as 

m? = ( -\cos2f3ml5 IJ + ml ij ), (4) 

with m e i (I = 1, 2, 3) being the mass of the charged lepton of the I-th generation. Two mixing matrices Z~ ^ 
which diagonalize the square-mass matrices of the sneutrino and charged slepton respectively are defined as: 

Z\m 2 Z- = diaqim 2 , m 2 , m 2 ) , 

Z[m 2 Z- = diag(m 2 , m 2 ,---,m 2 ) . (5) 

E E E E^ _Z?2 

As for the up- and down- type scalar quarks, we can define the mixing matrices similarly 

Z[ m 2 Z- = diaqim 2 , m 2 , m 2 ) , 

Z l m l Z D = dia 9(m 2 3i , m\ 2 , ■ ■ ■ , m\ g ) , (6) 

the expressions for those mass matrices m 2 , m 2 can be found in Ref.[30]. 

In order to suppress unexpectedly large effective FCNC interactions, it is natural to assume m 2 LlJ <C 

. Ax - . 2 

m 2 LlI , m 2 RlJ <C m 2 R u and A l u <C A l n with J ^ I. Accurate to order 0((^p—^-J ), we write down the 
expression of the mixing matrices which diagonalize the scalar fermion square-mass matrices 

Zj' J ~ U, (j , (i, J = 1,2,3) (7) 

and 

Zl ai = U. cos9. +U- sin#- , 

-E B oi E i E a(3+i) E i 

z t«(3+i) = _u. sinfl. e^i +U. cosfl. (a = 1, • ■ ■ ,6; i = 1,2,3) . (8) 

£ E ai E i E cc(3+i) E, 1 ' ' ^ W 

In general, the N x N transformation matrices U s can be written as 



|Ax s | 2 

1 S ij 1 



■ /■ 2(xq. x~,) 2 
Ax„ Arc, Ax, 

-y _ °ij _|_ ^ ' 'fc fcj ^g-j 

' J ^ ~ X S; fc^i, j ( X Sj ~ X S l )( X S j ~ X S k ) 

where IAto 2 I <C |m 2 — m 2 I (i, 7 = 1, • • • , AH. The symbols are defined as Ax„ = — t^-, x c = -V- 

with S = v, E. For the sneutrinos, Am 2 (z 7^ j) = m 2 ... The expressions for the off-diagonal elements of 
the charged slepton square-mass matrix are more complicated: 



Am 2 =cosfl. cos#- m 2 .. + sinfl. sin#. ^m 2 .. + 2m w g w c /3 



i<f> . . - 

+ 0036*- sin6>~ e , 



cos6»- sinfl- e 5 » A 1 *, 
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M3+j) 



Am 2 

E (3+i)j 



—i<f>B, o — i<f>~ 9 

— cos 6- sin 9- e im .. + cos6>~ sin#- e E im 



W-3 



Am 1 



Am 2 * 

B i(3+») 

n /i 2 •/!•/! "^F ) 2 

= cos y - cos 9 - m . . + sin - sin # - e 1 i m . 



4 (3+i)(3+j) 



2m w s w c /3 



cos0~ sin0~ e b j ' A^* + cos 6* . sin 6*- e B »A 



Am' 



Am 1 



i (3+i)i 



^i(3+i) 



0, (i, 3 = 1, 2, 3; i^j). 



(10) 



For a special case where there is degeneracy among the eigenvalues of the square-mass matrix, the explicit 
forms of the mixing matrices are given in appendix A. 

As we will find in following sections, the effective Lagrangian of e J e Ir y and e — > 3e are mediated by a 
combination of the following couplings 



G 
G 
G 
G 
G 



{ijigja} 
{ijigja} 

v m 

{ijigja} 
v (c) 

{i] 1 0. J a} 
v (d) 

{ijigja} 



+ + V V ' 

j^ii* 2, 2 i* 2,1^* Z^ a 
2^2i* 2^ ft* a 



m e i c w z (3+I)/3 z 3j 



C/3 



G 
G 
G 



{ijigja} 
Lib) 

{ijlf3.ja} 
£,(<=) 

{ijigja} 



2s w Zf^*Z- + (^w + ^) - ^^r^' 



^ V iv w ^ w 



m 7 c„ 



2-(3+J)a* ^3»" 2 S ^(3+-f)/3 jjlj* _|_ m e i C w gip ^3j* 



C/3 



2s ^(3+J)a*^H _|_ m e J C w yjq« 



w E N m w cp e " 



(11) 



Generally, we can recast the combined couplings in terms of Eq. (7) and Eq. (8). For example, we can write 

the coupling G , , as 

„(<0 



G 



{ijigja} 
,j(a) 



= z^*z l J 



Ax- i 

Sf3iS a j + 5 aJ 



VP 



Ax* 

+ S pr -" J 



x. — X- 



(12) 



where the first term just contributes to the muon MDM and EDM, while the other terms contribute to both 
the muon MDM, EDM and the FCNC processes of leptons. The other couplings can also be written down in 
a similar way, for saving space, we omit those concrete expressions here. 
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It is well known that one can apply the mass insertion approximation (MIA) to simplify the expressions 

of supersymmetric contributions to flavor changing neutral current (FCNC) processes which are induced via 

loop diagrams [31, 32]. In that approach, a small off-diagonal mass is inserted into the mass matrix which is 

written in the basis of the weak interaction and an approximate degeneracy of the squark masses is assumed. 

Thus the drawback is obvious while evaluating some processes where degeneracy of masses does not exist. 

Instead, in this work, we carry out all the calculations rigorously in the mass basis and keep appropriate 

pole masses in the propagators and mixing entries between various flavors at the vertices. We only expand 

the mixing matrix in the soft-supersymmetry breaking terms with respect to the mixing parameters 

/ Ax §-. \2 
2 _ l \ which are small as long as i f= j. 

\ X o. X a . J 

To be more explicitly, we would like to compare our approach with the MIA method, and point out the 
improvements of our scheme from the MIA. 

• When the mixing between left- and right-handed sfermions is negligible, i.e. S ~ (S = Ui, Di, Ei) 
and \Am 2 | <C \m 2 — m 2 \ (i, j = 1, • • • , N), our results are accord with the result of MIA at the 

lowest order of Ax s ... When S ^ (S = Ui, Di, Ei), indeed, our approach is an improvement from 
the MIA. 

• When all the eigenvalues of the mass matrix are approximately degenerate, just as proved in [32], our 
results are the same as that of MIA at the lowest order of Ax q . 

• As some eigenvalues of the mass matrix are only approximately degenerate, simple applications of the 
MIA method are invalid [33]. Now, there is a large flavor mixing in the sfermion sector (certainly, such 
flavor changing effects will be suppressed by the mass of the heavy sfermion in our detectable processes). 

In our numerical results, we will extensively discuss the first and third statements. In the following section, 
we give the formulae of the muon MDM, EDM, and the decay width of r — > fij. 



3 The muon anomalous magnetic and electric dipole moments 

The effective Lagrangian is extensively applied to evaluating rare decay widths of b, c-quarks[34, 35, 36, 37]. 
Derivation of the Lagrangian is carried out according to the principle: if all the masses m/s of the internal 
particles in the loops are much larger than the external momenta i.e. m 2 3> p 2 , thus the heavy particles can 
be integrated out. In our case, all the SUSY particles are integrated out and their contributions are attributed 
into the Wilson coefficients in the effective Lagrangian. 

For the W-boson propagator, we adopt the nonlinear gauge whose gauge fixing term is [38] 

gauge — fixing ^ (~^^) 

with / = (dfj,W +tJ- — ieA^W +fl — iS,m w (f) + ) and specifically we set £ = 1 in the later calculations. A thorough 
discussion about the gauge invariance in this situation has been given [39, 40]. 

The Feynman diagrams for e J e 1 ^ in the supersymmetric model are drawn in Fig. 1, the effective La- 
grangian is written as 

^ 7 = ^E^(/%)0?, (14) 
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and the operator basis consists of ten operators 



Of 



ot 



1 



(4vr) 2 
1 

W) 

1 



2 e J {itf,eF ■ a}uo T e ! , 



(4vr) 2 
1 



(4^2 , < 

1 



(4vr) 2 

1 



(47T) 



m je J eF ■ auj-e 1 , 
m eI e J ' {iFjlfu+e 1 , 
2 m I e J eF ■ auj + e I . 



(4vr) 2 ■■' 
1 



(15) 



This basis also exists in the case of SM [35]. In these operators, = da — ieQ e A^, F^ v = d^A v — dyA^ 
denoting the electromagnetic field strength tensor and F ■ a = F^ v <j^ v . The terms of dimension- four which 
are related to the e J ~f p uj±e I vertex cancel each other as long as we let e 1 and e J leptons be on their mass 

shells[35], so that they do not exist in £° e 7 at all. To shorten the text length, we present the one-loop 
contributions to the Wilson coefficients Cf{n w ) in the appendix. Here, we give the detailed expressions for 
the two-loop Barr-Zee-type diagrams (Fig.2). After integrating out the heavy degrees of freedom, we obtain 
the two- loop Wilson coefficients: 

etan (3 



a 



\n c £ Zf*Z^V KL (r H+ij ^y[Q u G{ XH _,x , Xi) ) 



+Q d G ( x H -> X n' X f 



y z Ki *Z K i*(T H+ ^Y G u x _ x _) 

Z-^l v E \ ) V H- ' Ej v i> 

K,i,j 



U,D 

x(T,+s 2 

\ i w 



T AQ*Qj( T y z Ki* z Kj _ Q 2 S 2 gij\ 



+ (t a q*AY(t yz Ki z K ^-Q 2 s 2 S'As 

/V'^QQ q w / 

+4V2yG(clx 2 A ,clx 2 ,c 2 x- )\T A ^^(T i yZ K »Z K ^-sl^ 

E R 

x (r, + 4) - (? A ^ j Y (r, E zf*T ~ s l 6ij ) « 

K 

+aV2 S ^\n c y QlT A ^F{^) + yT A ^^F{^)] \Su , 



U,D 



E 



c. 



i(Mw) = 



etan /3 



5(2-/oo P )VP W y i28A/2vr 2 m s 



if- ' Ui ' D; 
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K 



tf,»J 

W^V^c 2 * 2 ^ 2 ,c 2 x^ )[_r^(TV2^-QW) S 
U,D 

W2 £ G^x\ , cl*l , c^x ) [T^ ( Tf E Zf * zf _ s ^) s 
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X - 



X , 



>IJ 



(16) 



-4v& iV c J2 Q 2 q T AQ ^F{^) + Y J T AE ^F(- 
Xa u,d Xa e 

with q = m, di, Q Ui = §, Q di = \, T Ui = -T di = -T h = \ (i = 1,2,3) and N c = 3 is the color factor. The 
two-loop functions F(a), G(a,b,c) are given by 

'x(l — x) > 



F(a) 



1 , x(l-x) 
ax 



a — x(l — x) ^ a 

ax(l — x) In a 



G(a,b,c)=[ dxx{- — - — ^ 

Jo 4a-l)[ax(l 



+ 



x) — bx — c(l — x)] 
x(l — x)[bx + c(l — x)] 



In 



bx + c(l — x) 



} • (17) 



[ax(l — x) — bx — c(l — x)][x(l — x) — bx — c(l — x)] \ x(l — x) 

The couplings between the scalar quarks and Higgs are presented in the appendix C. Integrating out the 
heavy degrees of freedom in the loops, effective vertices of ~/HV where H is a neutral or charged Higgs boson 
and V stands for the vector gauge bosons W or Z, are obtained [25]. The gauge invariant forms of the vertices 
should be 

factor x [(k ■ q)g fiv - k^q v ] 

where and k u are the four momenta of the photon and vector gauge boson (W or Z) and the factor depends 
on the heavy degrees of freedom, the two-loop Barr-Zee type diagrams only induce the nonzero contributions 
to the coefficient Cf. 

Setting / = J = 2 in Eq.14, we obtain the muon anomalous magnetic dipole moment in the MSSM: 



SUSY GfITIu 

da u = i= „ 
V2tt 2 



ct + c 2 + l -(c^ + ct 



I=J=2 



Correspondingly, the muon electric dipole moment is 



,SUSY 

" ~ 6 i4vr 2 



I=J=2 



(18) 



(19) 



For the FCNC process r — > ^7, the contributions from SM and Higgs sector are suppressed by the small 

m 

ratio x v . = — j 1 . The supersymmetric contribution originates from the sneutrino-chargino loop and the 

* m w 
amplitude reads: 



(20) 
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where e is the polarization of the emitted photon. The form factors F L , F R are formulated as 



-F"k — (7o + Cc H —C<t , 



From Eq.20, we have the decay width as 

r. ... 



?R 



e 2 G%m 3 



m 

ci + c+ + ^C 2 . 
///. - 



128vr 5 



m 



+ m 



?R 



(21) 



(22) 



4 The lepton-flavor-violating decay r — > 3/i 

The effective Lagrangian for r — ► 3/u. is induced by the following four pieces: 7, Z, .ff-mediating penguin and 
box diagrams. Those Feynman diagrams are shown in Fig. 2 and Fig. 3. After integrating out these heavy 
degrees of freedom, the effective Lagrangian is written as 



T— >3/U 

eff 



G%m 2 x 



(23) 



with those operators being 



Q4 = jiU-TjlUJ+H , 
Qq = jMjJ+TpUjJ+H . 



The differential width of r — > 3/x is 



d 2 r 



dm^dm^ (2vr) 7 32m,3 



l-MI 2 , 



(24) 



(25) 



where mfj = (pi + Pj) 2 , and pi (i = 1, 2, 3) are the momenta of the outgoing muons in the rest frame of r. 
With Eq. (23), one obtains the square of the transition matrix element \M\ 2 as 

\M\ 2 = {4(m 2 + m 2 - m\ 2 ) (m\ 2 - 2m 2 ) (|d| 2 + |C 2 | 2 ) 

+4(m 2 + m 2 - m 2 3 ) (m 2 3 - 2m 2 ) ( | C 3 1 2 + |C 4 | 2 + |C 5 | 2 + |C 6 | 2 ) 

+8m T mjRe(4CiC 2 t + C 3 Cj + C 4 Cj) + 4m 2 (m 2 + m 2 - m 2 2 )Re(c7iC] + C 2 C\) 

+4™* (m 2 + m 2 - m 2 3 )Re(c , 1 C 4 t + C 2 C\) 
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+4m T m M (m 2 + m\ - m\ 2 - mi 3 )Re(CiC 5 t + C 2 C]) 
+4m T m^(m 2 12 - 2m 2 i )Re(c 1 C 6 t + C 2 C\) 
+4m 2 (m 2 + - m 2 3 )Re(c 3 Ci + C 5 C\) 

+4m T m^ (m| - 2m 2 )Re(c , 3 C , 5 t + C 4 C{) } . (26) 

5 Numerical result and discussion 

In this section, we present our numerical analysis on the muon MDM, EDM and the lepton-flavor-violating 
decay processes in the supersymmetric scenario with the nonuniversal soft-supersymmetry breaking. In the 
lepton sector of MSSM, there are 15 + 3G (G denotes the generation number) new free parameters besides the 
SM parameters gi, g 2 , m 7 (I = 1, 2, 3). Too many parameters reduce the predictability of the model and 
before any direct evidence of the SUSY particles is found, determining the free parameters is the most subtle 
and tough task in the supersymmetric theory. To find a way out, the grand unified theory (GUT) assumption 
is frequently adopted where all new physics parameters are fully fixed from only five free parameters at the 
Grand Unification scale. In that scenario, all the input masses conserve flavors, namely there are no off- 
diagonal masses are to be introduced. Thus when all the concerned parameters evolve from the GUT scale 
down to the electroweak scale by the renormalization group equations (RGE), the flavor changing parts in 
the effective Lagrangian emerge only through the CKM mechanism. The induced FCNC-related decay modes 
would have very small branching ratios which are much lower than the experimental bounds. Therefore, we 
would rather adopt an alternative parameterization which may imply a different physics picture from the GUT 
scenario, i.e. to choose a set of parameters at the electroweak scale as inputs. The numerical results depend 
on these parameters and by imposing the experimental bounds, one can obtain constraints on the parameters. 
Just as in Ref.[27], we require the phase arg(n) = in order to suppress the one-loop contribution to the 
electron EDM. At the lowest order we can choose the parameter basis which is responsible for changing lepton 
flavors as following: 

H, mi, m 2 , tan/?, m 2 ^, rn 2 _ ^ > ^ ^ m \u' A \j' ( J > 3 = lj 2 ' 3; 1 ^ J ^ 

Correspondingly, the square mass of sneutrino is given through the relation: 

m 2 = cos 2 6 - m 2 + sin 2 6 ~ m 2 — m 2 + m 2 cos 23. 

V I E I Ej Ej e (3+i) e i w 

There are several relevant CP phases as : 4>/i = (frgi for I = 2 4>f = (fi^i for 1 = 3, and 4>a\^ 
m 2 , § m 2 R {I 7^ J) are the CP phases of parameters Ajj, m 2 Lij and m 2 Rij respectively. To simplify 
our discussion, we assume that the bilinear and trilinear couplings of scalar quarks are all universal, i.e. 
Af T = A rSjj, A d Ti = A ,5u and m 2 = m 2 5jj, m 2 = m 2 <5/j, m 2 = m 2 Sjj. As aforementioned, 

the nontrivial CP phases lead to a large mixing among the neutral Higgs fields [29]. Considering the two- loop 
Yukawa and QCD corrections to the effective potential, the square mass matrix for the neutral Higgs is written 
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as: 



m 



ho 



2 2 



2 2 
W W 



+ Re(A 5 )s^ +Re(A 6 )s (3 c^ 



/ 



-m„ s , 3 c , 



a 2 s 2 
W W 



°/3-/3 J2" 
+ Re(A 6 )c 2 + Re(A 7 )s 

\ Im(A 5 ) S/J + Im(A 6 )cg 



A 3 + A 4 



a 2 s 2 
w w 



+Re(A 6 )c 2 + Re(A 7 )s 



/ 



„ 2 2 

m 2 c 2 W_W 

a 



+ Re(A 5 )c 2 + Re(A 7 ) S(jC(3 



/ Im(A 5 ) S(3 \ 
( +Im(A 6 ) C(3 J 



/ Im(A 5 ) S(3 \ 
^ +Im(A 6 ) C(3 J 



Im(A 5 )c fl + Im(A 7 )s £ 



(27) 



with the square mass m 2 



2 2 

m„ = m , 



4m, 2 „s 2 „ 



^A 4 - Re(A 5 )) . 



,- - / < 28 > 

Here, the parameter m ff± represents the mass of the physical charged Higgs-bosons, and the concrete ex- 
pressions of the other parameters Aj (i = 1,2, ■■■,7) are presented in appendix F. As pointed out by the 
authors of Ref. [29], the experimental lower bound on the mass of the lightest neutral Higgs mass can be 
reduced to 60 GeV in the CP-violating MSSM. In the numerical analysis of this work, we take this lower 
bound as an input. With above specification about the parameter space, we carry out our numerical compu- 
tations. Without losing generality, we take m R± = 300 GeV, m\ = 1 TeV, m gi = m So = 10 TeV, m fii = 



m fi = 500 GeV, m- 2 = m. = 700 GeV, 6 S = <f> e = 0, 0- = 6 f = \, = m\ w = = 



GeV 2 , A[ 2 = A[ 3 = A l 21 = AL = GeV, M s 



m 



l 12 

500 GeV 



13 

A t 



1 TeV, m t = 174 GeV, m b = 4.5 GeV, 



rn 



Q 3 



m-„ 
u 3 



e*4 TeV all through the paper. Indeed, the electron EDM is a very rigorous 



21 

f,3 — — . , — Af, 

constraint to the parameters. Although we set the CP phase which is only related to the first generation of 
sleptons to be zero (in this case the one- loop contribution to the electron EDM is also zero), the non-zero CP 
phases which are related to the second and third generations can also lead to a large EDM of electron via 
the two-loop Braa-Zee diagrams. For some regions in the parameter space, the two-loop Barr-Zee diagrams 
would contribute an EDM of electron which is larger than the experimental upper bound. Thus the electron 
EDM restricts the relevant parameters via the two-loop diagrams. Considering the constraint of the electron 
EDM through two-loop Barr-Zee diagrams involving sleptons, we set the CP phases <f>- = (f>- = 

At present, the experimental upper bounds on the branching ratios of the two lepton-navor-violating 
processes are 

B(t -► M7) < 1-1 x B(t -► 3/i) < 1.9 x 10~ 6 . 

Upon imposition of the experimental limits on our theoretical predictions for the branching ratios and the 
muon MDM, we find that the contribution to the muon MDM from slepton generation mixing parameters is 
much less than the contribution from the flavor conserving parameters through scanning the parameter space. 



Taking 



i2 
«23 



and \m 2 r 



l L 2 z\ ~ l m /?23l — ±uu ^ cv ' 1^231 — 1^231 

plot the muon MDM versus the parameter TO2 in Fig. 4. Within la tolerance, we find that the theoretical 
prediction coincides with the experimental data rather well if a suitable parameter range is adopted. As for the 
calculation on the muon EDM, we must consider the constraint from the electron EDM. With 2a tolerance, 
the experimental upper bound on the electron EDM is \d e \ < 0.5 x 10~ 26 e • cm. With mi = 500 GeV, 4> , = 



100 GeV 2 , \A l 23 \ = \A l 23 \ = 100 GeV, we 



A 1 
A 32 



^23 



2 ' 



i L 3 l 



\m 



R23 1 



100 GeV 2 , |^ z 23 | = \ A l 23 \ = 100 GeV, and tan/? = 5, 50 



respectively, we plot the electron EDM, the muon EDM, and the muon MDM versus the parameter fi in Fig. 
5. We find that the muon EDM can reach 10 -23 (e • cm), which is one order above the proposed sensitivity of 
the coming experiments [41]: 10 -24 (e • cm), while the muon MDM and the electron EDM are also consistent 
with the experimental data. 
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Now, we analyze the lepton-flavor- violating decays of heavy- lepton r: r — > fij, r — > 3/x. With /z = 200 GeV 
and set all CP phases to be zero, we plot the branching ratio of r — ► 3/z versus = = 5m 2 (Fig. 

7(a)) and ^4 2 3 = ^23 = ^ (Fig- 7(b)). This result is slightly different from that given in previous literature. 
We not only consider the contributions from the 7— penguin, Z— penguin and the box diagrams, but also 
include the Higgs-penguin. We notice that for larger tan f3 values, the contribution of the Higgs-penguin is 
not negligible, It is in analog to its role for the rare leptonic decays of B-meson. Numerically, we find that the 
contribution of the 7— penguin is much less than that of the Z— penguin and box diagrams. For the case of 
non-zero A l 23 ,A l 32 and large tan/3, the Higgs-penguin is also non-negligible. When A l 23 = A l 32 = 0, but m 2 L23 , 

and fn 2 are not zero, the main contribution to the decay width of r — > 3/i comes from the the Z— penguin 
and box diagrams. By this we can explain why the difference of the curves corresponding to tan f3 = 5 and 50, 
is rather small, whereas in Fig. 7(b) the difference is so obvious. As to r — ► fij, when we take into account the 
constraint from r — > 3/i, we find that its branching ratio is obviously lower than the experimental bound. If 
the lepton flavor violation originates from the term m 2 , m 2 (Fig. 7(a)), the difference between tan/? = 5 

and tan/? = 50 is very small. If the lepton flavor violation originates from the term A\j (Fig. 7(b)), the 
difference between tan (5 = 5 and tan (5 = 50 is about two orders. We also investigate the r — > ^7 in the 
model. Imposing the constraint on the parameters from r — ► 3//, we find that the branching ratio for r — > ji r ) 
is smaller than 10 -8 , beyond the present experimental detecting ability. Assuming that the flavor violation 
originates from m 2 .., m 2 .., we plot the B{r — > fi'y) within the possible parameter range in Fig. 8. The 
situation is similar to the case of r — > 3//, namely if the flavor violation originates from the term A[j, the 
branching ratio B(t — ► ^7) is much less than that if the mixing m 2 ^. , m 2 induces the flavor violation. 



6 Conclusion 

In this work, we investigate the muon anomalous magnetic dipole moment, muon electric dipole moment, and 
the branching ratios of r — > fij and r — > 3/U in the framework of the CP-violating MSSM with the nonuniversal 
soft-supersymmetry breaking. For the muon anomalous magnetic dipole moment and electric dipole moment, 
the main contribution comes from the parameters which conserve the flavors. Process r — ► 3/i occurs mainly 
through the Z-, Hf—, A - penguins and box diagrams. It can help understanding why the branching ratio 
of r — > /X7 is so small as long as we consider the constraints from r — > 3fi. From the methodology aspect, our 
method is equivalent to the MIA scheme when all mass eigenvalues of the mass matrix are almost degenerate, 
or the square mass difference among the eigenvalues are much larger than the corresponding flavor violation 
parameters. For the case where only several mass eigenvalues are degenerate, the MIA is invalid. Indeed, our 
method is an improvement of the MIA. 
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A The mixing matrix for the sfermion sector 

For the general case, the mixing matrices for the sfermions are given in the text. In this appendix, we present 
the mixing matrix for the case where there is degeneracy in the mass spectra. In fact, we can always perform 
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a transformation on the sfermion mass matrix and turn it into a standard form: 



Z„ m a 

b LR S 



2 Z 



1 

S LR 



m 2 Si 
Am 2 * 
Am 2 * 

s 13 



Am 2 * 

b l5 

Am 2 * 

S 16 



Am 2 

2 

S 2 

Am 2 * 

£>23 

Am 2 * 

b 2A 



Am 2 * 

s 26 



Am 2 
Am 2 

s 23 


s 3 

Am 2 * 

6 34 

Am 2 * 

b 35 







Am 2 

Am 2 

S34 

m s 4 
Am 2 * 

6 45 

Am 2 * 

s 46 



Am 2 

s 15 



Am 2 
Am 2 
m% 
Am 2 * 

s 56 



with Am 2 being given in Eq. (10) and the matrix Z s is formulated as 



Am 2 
Am 2 

s 26 



Am 2 

*46 

Am 2 

J 56 
s 6 



m 



( 


cos 9„ 








sin#„ e l( ^ s i 










cos 9„ 

b 2 








sin^„ e 










cos 9 q 

S 3 










— sin 9„ e^ s i 








cos6L 










— sin#„ e s 2 

s 2 








cos 9„ 

s 2 


\ 








— sin#„ e^ s 2 

S 3 









sm( 






, e 

'3 




cos# 



Without losing generality, we assume that the first n (2 < n < 5) eigenvalues satisfy the condition m 2 — 
Am 2 (i, j = 1, • • • , n), for the scalar lepton sector the mixing matrix is given as 



with 



and 



Z s = Z s Z.V 

b b LR A 



|Ax„ Ax, I 2 \Ax 

Su ~ ~ 2—i jz _ 1 \2jz _ rf \2 ~~ 2—, tz z -T- "12 



Ax c Ax„ 1 |Axo | 2 Ax Q Ax 

-y _ \ ^ b ja b ai ^ \ ^ 1 b ai 1 b jf3 b 0i 

Sji ~ 2—, („ _„ \(„ _ ^ ^ („ _ „ \2 2—! 



a ( X Si X S j )( X S l X S a ) ( X Si X Sj) 2 a ,/3 ( Xs i X S a )( X Sg 3 



Ax„ Ax, Ax Q Ax,, Ax„ Ax,, 

Y _ s ai _|_ y^ s <*p s pj |_ y^ y^ s «j s jp b fa 

Sai nr nr ' J f nr v \( nr nn \ / J / J 



X Si X S a f3^ a ( X Si X S a )( X Si X S p ) j p ( X Si X S j )( X S i X S a )( X Si x s p ) 

Ax„ Ax„ Axo 



x s a x Si a^ n ( x s a x Si)( x s a x s„) 



v s = 1 - Y, 

13 



( x s a x s„y 



Ax Q Ax„ Ax„ Ax,, Ax, 



Y |_ y^ a /3i b ia |_ y^ 



( X ci X o ) ( X o X a ) ( X q X o ) / n ( X <~i X o X a ) 



I A 

I nxn nx(6-n) 

1 1 

\ (6-n)xn (6-n)x(6-n) 
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Here, the matrix A is used to diagonalize the block matrix m 2 (i, j = 1, ■ ■ ■ , n). Expanding the Eq. (31), 

ij 

we have 



where 



Z\ ai = U . cos 9 - + U . sin 6 - e , 

z t«(3+i) = _u sin # e ^ + u cos e ( a = 1, . . . , 6; i = 1, 2, 3) , (34) 

k 

u s . =v s . , 



U s =VV S A s , 

o a i t—^< a ak a ki 

k 

U, = V„ . (35) 



B Expressions of the coefficients in Eq.(14) for the lepton anomalous 
magnetic dipole moment and the amplitude of decay mode e 1 —> e Jr f 

The expressions of the coefficients in Eq.(14) for the vertex e I e J 'j are 

cr(^w) = - 2eG{ X" Jalx l- B Z]( x K 7 ,x * a " > 

i i 

1 {iilaJa} 

w 1 



1 {ii/aJa} 



4c 2 i (a) 
w 

2 {ii Ja J a } 



3 » 

7 „2 »«» 



(o) - \ 5 „(o) 



pi] v "°' 2 «" I 31 ' «■ 

1 {ii/aJa} (2) 



g c 2 L (a) [41] V b«' K 



" W 



C 4 (A*w) = -^^r^^r^^iK-'^) 

1 ^{iilaJa} i (0) 

w 1 

w 
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1 {iilaJa} , . 

w 1 

1 {iilaJa} . . 

^i(» fi a»» o). 



4c 2 i< d > 
w 

C 3 + (Mw ) = " ^2 ( X .r * e J ) * K (*,- . ^ ) " ^ (\- , ) 



7 .2 „(0) 



1 {ii/oJo} (2) 



% % W 



V2. 



W ' 1 

C 5 + (m w ) = -^G i 2; Ja \x K _x^B^ ] {x K _,x oa ) 

+ ^ G T(o) Ja}X E^ B Z ( ^ ' X ^ ■ (36) 

w 

C The couplings among the sfermions and Higgs 

In this appendix, we give the couplings among the scalar fermions and Higgs fields that involve the two-loop 
Barr-Zee-type contributions to the Wilson coefficients: 

2 2 

m , m 



^ "V 7O+ I)i 7(3+ J).?l + / e / i * m „J ylJ* _ yKJ* A u c \ Z {3+I)i z Jj 



+ — 

m 



\^/2rn « KJ 0) v D J 



m,„s 



2 

em 



Y^Uiiij = A em ui f^z H *Z (3+I)j - fi*Z Ij Z^ +I ^*) + — (A u Z Ij ^( 3+J )" _ z u* z (A+j)j\\ ^ 



pADi^j = A fj_( A d* z Ij z (Z+.J)i* _ A d z Ii* z (3+J)j\ _ em d i ( gli* Z (3+I)j _ * z Ij Z (3+I)i*\} 
I \ IJ D D i.i i i j 2m w s w \ D D ■ /J 
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T AE t Ej = A b P Z Ij Z (3+J)i* _ A l z Ii* z (3+J)j\ _ e "V (uZ Ii *Z (3+I]i - u*Z /j Z (3+/)i *H 

ly^V IJ EE I. J E E ) 2m... S... V EE r 1 E E J J 



2m w s w 



D Expressions of the coefficients in Eq.(23) for decay mode e — > 3e 

The contributions from the gamma penguin diagrams are written as 



4s 



W q(«'"M 



3 *<») 



, ( x ■ x ~ ru i x 13 , , (x . X' ~ 

pi] v 2 k~ ' 31 1 «~ " 



-^* 2 S (0) ,(x , x . a )l+iH_G { " J ; Ja} iBf 1 (x. 
2 K - I«P «r' " ay J 3 C 2 L (o) [4i] v e K o, 



-I 2 «(°) / 

X IS) . . I X « tZ» ~ r 

2 «: «r' - c 



3 C 2 

w 



B, Ax - , x „ 

[41] ^ £«' k9 



C£(// w ) = -2C?(// W ) , 



C 3 7 (Mw) = C 6 7 (/u w ) = • 



The pieces from the Z— penguin diagrams are 



cos 2 2(9 



^2(Xn,X- 



COS 26> w {ij/aJa} 



4c 2 



+ ^ X K X K ) 2 ) X K > X K , 



f z Ai* z Aj _ Z Si Z 3j*\ 
\ N N N N J 



It (1) I \ 

9 [in] ' x K o ) x K q J 



COS 2fl w {ii/jgja} ^ 2 

2c 2 



w 

2 COS 2 26 w {iilaJa} 

n ^ 



+ cos20 w G 



(a) 



f 2s 2 ^ - z^Z^*) T (x - ,x. B ,x ) 

V W EE) [HI] V -E Q E' 3 kV ' 

i 

-T W (x. a ,x ,x )(25 ij cos26 w + Z li Z lj *) 



-(x.._s _)27[ 111] (s. a ,x (( _,s (( 



2<^'cos20 w + 



+ cos20 w G , s Z (x oa ,x„,x _) , 



C 2 
w 
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\ N N N N ) 



-( \W (0) ( 

V X X 0/ i [ 11 ll \ X Ra > X K ' X K 



2c 2 

w 

K I [111] V E a 

j 

r 2 \ W EE/ 



{ij I a J a} 
L(d) 



-T W ( \ 

2 [mi ' x K o ' X K ° I 



{UI0Ja} (1) , . 
L (d) i [lll]V X B"' X B/3' X «oi 



-^r (* / x fi j ) * G ^7 ^ ( x .- ' x ^ ) 



--f(x e/ x eJ )fG 



1_ {jj/a Ja} 



It (1) c 

2 i [in] ' X K - ' x „ 



2<5^'cos20 w + 

- {x^ ) §7^ (x, a , X k _ , ) (2^" cos 2# w + ZlVZ^ 



2s 



C 2 



\ x „i x r j) 2 [in] \ x a a ' x ,-,8 ' X „ - > ' 



Cf(Mv 



Cf(/Av) 



As 2 
cos 20 z 



C#(/i w )=CfGu w )=0. 

The CP-even Higgs penguin-diagram contributions are formulated as 



<?f fc Vw) 
^(Mw) 



Cf fc °(/x w )=0, 



1 



2c 2 c 2 ft x n 



{iila Ja} 



+ 



J [ni] y x E<* > x K o ' x K o ) \* R ^ N ^ R ^ N 



1 \ 7 ifcp{^ /Qja > 



(^ S w-4\)+Kf/C]( 



« J 



T 2 (x K7 ,x oa )(z™) 2 [(x eI x^G { 2 



1 {ii/a Ja} 



+ 2 C 2 ,<««) 



(39) 
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2 . .1 {jilaja} 

i X r J X r ,j) 2G 



2 

c /3 x H o 

k 



+Z 2 R k Z 1 _ i Z 2 + j ) 



li* 
+ 



) + < 



2i 7 lj 



rlk 



2^2' 



m w" s w C w C I X jfO 

k 



-(x^xJ^B^^B^ix^xJ, 



cf*Vw)=^ fc Vw), 



1— ^(x^.x-J^) 5 



{iila Jot} 



_L {iiJct Jq} 



i 2 _s jri/tr 1 ^™ -7 "* 



+ ( x I x j) 2g , , 

T [lll] ( X E<* ' X K ' X K ) ( Z R kZ N ~ Z2 R Z N 



(*x-^)+(*^o)*0- 

» j ■ 



x- ,x „,x „ll£^ 



k 



«( v ,x^)(^)*[^G^-" + 



9 ^ I 

C P X h° 

k 



T m (x. a ,x ,x )(z lk Z 2i Z 



+ 



+z 2k z 1 J*zl i * 

R T" 



2-\/2m„ .s„ r- r-.r 



,2 „2 



x ^_G iU ': Ja} Z lk B k B { "\( x - a , x ), 

e J k . „fcl R R 21 v " ' 



Ce fc (Mw) = C , 5 (M W )- 



(40) 



The CP-odd Higgs contributions are 

cf(» w ) = cf(» w ) = o, 



C£ (Mw) = 



2c 2 x ^ tan ^ 2 (^„,^)[(^X e ^G (i) +^ G i(a ) 
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2c w x aO 



l—xl tan p [ - (x, Q , X k0 , X k0 ) (spZ? - cpZf) 



+— 



G 



{ijlaja} 



H tan /?JF 2 (x _,x. a ) 



{ii/aja} e I X e J ^{iilotjoc} 

X < J »W 2cl 



(x eI x £j ) 2 tan 2 /3 (x.„ , x^ , x^_ ) (spZ^Z. 



s P x A o 



2j* 



+c (3 z 1 J*z**) - (x k7 x k _ ) ^ ( Xpa , x k7 x kt ) {sezpz? + c^zf) 



' j. 2 a<r- I \\ ~,{UIaJa} . ,l p {aoJ«) 

k tan pJ-oyx n ,x~ ) x T G , , + (x r x . 2 

2 C 2 x ^ L ^ L d ^ e' e-J> L (a) 



G 



{jilot Jot} 



"W X AO 
1 



4 tan /3G^ Ja} (x^ , , X r0 ) ( 8fl Z* - c^j) 



(4' s w " 4' c w) - ( I „o^o)^p° 1 > 11 (^.^>^)(^ 



3i* 

JV 



JV 



J-tan 2 /3^ 2 (x ,x„j[(x ,* ^G^ + ^^G^ 



AO 



h8 



^ao 



-x . tan pG , s 

„(c) 



-(x^x^_)2^ iii] (x 5a ,x^_,x 



The box diagram contributions to the coefficients are 

r box ( {hipjo} {jj j a m (i) 



w 



{ii//3J^} {jjJotJa} (1) , . 



(41) 
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c box ( {mpj<*} wjcjpy CD , , 

°2 (Pw ) ~ 4( ,4 ^ i(d) ^ (<J) ^[lll] l X K > X K . X £a . ^ 

1 {UI0J0} {jjJcxJa} i i (1) 

W 1 3 



2 {UI0J0} ^{jjJceJa} , N 1 ^(0) 

-/3 



C box { J_r G {«//»J-} G {«Ja^} 1(0) , x 

W 1 3 % 3 

1 {<i//3Ja} {jjJaJ0} (1) "I 

2 z» z,(<9 l 111 l v «° «?' B " b /3; J 

— 1 J 

, 2 (~,{jU0J0} {jiJaJa} . ,1(0) , . 

H TiCr , Lr , X ,(X X ) 2 L> , AX ,X , X - a , X a ) 

°/3 1 3 1 3 

{iilpjp} „ {jjJocJc} X j (1) 1 

— Ljt Cj , x -^—L), Ax ,x ,x. a ,x a )> , 



r 6ox, N _J_r G WJaJ-} {jiJ/W}, 1(0) 

w 

1 {iUPJcx} {jjJaJ0} (1) "I 

2 i(d) z» [HIP «?' E°' S/3 ^ J 

, 2 f < ^,W//3^}„{ijVaJa} ,1(0) , , 

+ ■%{ G (*) G ( d ) *V X K - X K ~ ) 2 P [iii] ' \t ' x *» ' X ^ ) 



X E a ' X B/3' 



~{jiW/3} „{iW«Ja} >, X e / X e J^ _(i) , O 

+ d G ,W G ^) ) 2 ^[111] (^ K r ' X .7 ' X - ' X ^ ) 

^(3 1-3 13 



(42) 



E The loop integral functions 

In this appendix, we present the expressions of various loop integral functions which appear in the text. 

x\x%Vca.x\— lnx2) x 2 — 5xiX2— "2xi , £ , , 

(ss-si)* 6^ '( fOT ^1^^), 

i^i . ( for »i = ^2); 



/i(a;i,X2) 
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f2(x 1 ,X 2 ) 



B™ ] (x 1 ,x 2 ) 



R (0) ( \ 



BZ(x 1 ,x 2 )=BZ(x 1 ,x 2 ) + x 1 BZ(x 1 ,x 2 ), 



(2xiX2— x|) lnx 2 — x 2 lnxi 
2(x 2 -xi) 2 



El 



gi 



2(x 2 — Xl) 2(x 2 — Xl] 



, (for X 1 / X 2 ), 



•j - ^lnxi , (for xi = x 2 ); 



i 

' 2xi 



(for xi = x 2 ); 



6x 



^ , (for ari = x 2 ); 



= < 



2(lnxi—lri3:2) _|_ 2x 2 +5xix 2 



(x 2 -xi) 4 



6xf(x2— xi) 3 



, (for xi / x 2 ), 



12^ , (for xi = x 2 ); 



,(0) 



(0) 



[31] 
-,(2) 



[21] 
,(0) 



[31] 



B)"' (xi, x 2 ) = BT. (xi, x 2 ) + 2xi5 ( ° ) l (xi, x 2 ) + xifB'"' (a?i, z 2 ) 



2 D (°) 



[41] 



^(xi, 0:2,2:3) = < 



[21] 



[31] 



[41] 



xi In xi 



X2 lnx 2 



X3 In X3 



(X2-X1XX3-X!) " (xi-x"2)(x 3 -x 2 ) - (xi-x 3 )(x2-X3) ' ( fOT X l ^ X 2 + X 3 ), 

X3 %:-xP ] + ^7 , (for x 1= x 2 /x 3 ), 



X3-X1 



> -2iT > ( for x i = x 2 = x 3 ); 



^| 1 1 ) 1] (xi,x 2 ,x 3 ) = { 



Xjlnxi x?,lnx 2 x|lnx3 /,, _L _L \ 

~(x 2 -xi)(x 3 -xi) ~ (xi-x 2 )(x 3 -x 2 ) ~ (xi-x 3 )(x 2 -x 3 ) ' i tOT X l ^ X 2 r X 3 J, 

(2^gnxi^ + _^ i(for Xl=X29 , X3)j 



X3-X1 



P [llll] ( X 1' X 2,X 3 ,X4) 



lnxi , (for xi = x 2 = x 3 ); 

_ xi lnxi x 2 In X2 

(X2-Xl)(x 3 -Xl)(x 4 -Xl) (X1-X2)(X3-X2)(X 4 -X 2 ) 



X3 I1IX3 



X4 In X4 



(xi-X3)(x2— X3)(X4— X3) (xi— X4)(X2— X4)(X3— X4) 



, (for Xi / X 2 / X 3 / X 4 ), 



X3 In X3 1 X4 In x 4 

(X1-X3) 2 (X4-X3) (X1-X4) 2 (X3-X4) 

~(x 3 -xi) 2 (x4-xi) 2 ~~ (x 3 -xi)(x4-xi) ' ( for Xl = X2 ^ X 3 / ^4), 

_2xi]n^^lnx a _ ? (for ^ = ^ ^ Xg = ^ 



(x 3 -xi) 3 

X4(lnxi— lnx 4 ) 
(X4-X1) 3 



+ 



(x 3 -xi) 2 

2 ,. i ,,. 1 _ ri) . ,(for xi = x 2 = x 3 /x 4 ), 



X4 



^ , (for xi = x 2 = x 3 = x 4 ); 
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x\ lnxi 
(x2— x\){xz—x\){xA—x\) 
x\ In 23 

(x\ -X3) (X2 —X-i) (X4~ xz) (x\—X4,)(x2~X4,)(xz-X4,) 



x\ \l\X2 
(xi-X2){xz-X2){x4-X2) 



x\ In X4 



, (for X 1 / X 2 / X 3 / X 4 ), 



x|lnx3 x\\\\X4 
(xi-x :i )' z (x4,-x 3 ) (Xl-X4) 2 (x :j -X4,) 

(^xxxzx^-xixz-x^x^Xnxx x , , f _l _l \ 
{x3 - Xl ] Hx4 XY {xz-x^U-x,) > ( fOT «1 = ^ / X 3 / X 4 ), 

2^jnxi-2x|ln^3 +23+2x1 In ^+2x3 lnavj) , f _ / _ n. 



(43) 



+ 2(l"V > ( fOT ^1 = ^2 = X 3 / X 4 ), 



x|(lnzi— lnz4) 



(a; 4 -a;i) 3 

-3^ , (for Xi= X 2 =X 3 = X 4 ). 



F The expression for the parameters 

The A parameters that involve the neutral Higgs mass mixing are given as 

x . 91 +92 3 2 X 



16vr 2 



hi 



1 + \ Xh + lf^(l^ + ^ " 85 {Xbt + 



+ 



i92tt 2 *M S 4 USY 



1 + 



16vr 2 V 2 
ii (9/l 



16 5s 2 )t 



16vT 2 ^ 



+ 



* + \* + + \hl ~ 8g 2 ) » + t 2 ) 

M 4 



192vr 2 b M| USY 



1 + 



16vr 2 



(9h 



5*? 



16<tf)t 



A, = - 



52-5i 



16vr 2 



r 1 




l t+ 2 


Xtb + 1 6 7T 2 


3|H 2 


M 2 I^| 2 


MflXSY 


M SVSY 


3M 2 


M 2 I^I 2 


M SVSY 


^SUSY 



1 + 



1 + 



1 



16vr 2 
1 

16^ 



(6h 2 - 2h 2 - 15g 2 s )t 
(6h 2 b - 2h 2 - 16g 2 )t 



9l 



1 



3 

ST 2 



16vr 2 

2l2 



(h 2 +h 2 )t 



+ — W 



1 1 

t+ 2 X * 6+ 16^ 



(ft? + fcg - 8<? 2 ) (AT t6 i + t 2 ) 
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h 



96vr 2 * V Ml 



SUSY 



96vr 2 



3|/[/| 



M SVSY 



|^| 2 |^| 2 
/If 4 

JW SUSY 

\»\ 2 \A b \ 2 

^SUSY 



1 + 



1 + 



1 



16vr 2 
1 

16vr 2 



(6h 2 



2h 2 - \Qg 2 s )t 
(Qh 2 - 2h 2 - l&g 2 s )t 



3 h A £A 
192vr 2 *M| USY 



16vr 2 



( 2/l 2 _ 6h 2 t + 16 5s 2 )t 



+ 



1927T 2 



a 2 A 2 

iw SUSY 

3 4 \iA 2 yA t 

96vr 2 * M 4 USY 



1 



16vr 2 
1 



(2/* 2 - 6/i 2 + lQg 2 )t 



16vr 2 V 2 



(\h 2 b - ^-h 2 t + 16g 2 s )t 



+ -A- /j 4 ^ f _ 

96vr 2 6 Msusy V Msusy M| 



A| 2 A 



A 7 = - 



+ 



96vr 2 
3 



^SUSY 



1 



16vr 2 



SUSY 
2 * 



6A 



\A t \ 2 A t 



96tt 2 ' * Msusy v M SU sy 
where t = ln(Mf USY /m 2 ), gi = -£-, 92 = 3^ and 

V2m t (fn t ) 



^SUSY 



/ L 16vr 2 V 2 

- Y^6 + i% 2 )^ 



/i 2 + 16 ff 2 ) 



^ 2 + 16fl2) 



/it = 



v sin/3 

2|A f | 2 
M 2 

JW SUSY 

2\A b \ 2 
M 2 

iw SUSY 



h b = 



\A t \ 2 
12Mg USY 

\A b \ 2 
12M| USY 



V2 



m b {m t 



v cos j3 



)• 
)• 



|^| 2 + |A b | 2 + 2Re(A*A f ) 

9 A/f 2 

z iw SUSY 



I /i 



: 2 - AjA t 1 2 



M 2 

JW SUSY 



6M 4 USY 



In Eq. (45), rat is the top-quark pole mass, which is related to the on-shell running mass rrit throug 



m t (m t ) 
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Figure 1: The Feynman diagrams which contribute to the e e 7 effective Lagrangian in MSSM 
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Figure 2: The two-loop Barr-Zee-type Feynman diagrams which contribute to the e ' e '7 effective Lagrangian 
in MSSM 
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Figure 3: The Feynman diagrams which contribute to the e — > 3e in MSSM 
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Z, Hi A \ Z, Hi A 




Figure 4: The penguin diagrams in MSSM where Z ,H®, A are involved respectively 
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Figure 5: The muon anomalous magnetic dipole moment versus m.2 in MSSM with (a)/i = —200 GeV, 
(b)/i = 200 GeV and tan /3 = 5, 50, the other parameters are set as in the text. The shaded region is allowed 
by la tolerance from the most recent experimental observation. 
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Figure 6: The muon's MDM and EDM, the electron's EDM verus the parameter \i in MSSM with (a)tan (3 = 5, 
(b)tan (3 = 50; the other parameters are the same as in the text. 
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Figure 7: With \x = 200 GeV and all zero CP phases, the branching ratio of r — ► 3fi versus (a) m? L 
m «23 = or ^) ^23 = ^23 = the other parameters are taken as in the text. 
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Figure 8: The the BR(t — > a-y) versus m? = m? = 5m 2 , the other parameters are taken as in the text. 

v r " L 23 fi 2 3 



34 



